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1.   INTRODUCTION 

The  year  1986  is  likely  to  be  recorded  in  the  world  history  as 
the  time  of  the  greatest  nuclear  disaster.  With  the  public  aware 
that  oil  supplies  are  limited  and  nuclear  accidents,  like  the  one  in 
Chernobyle,  yet  possible,  increasing  emphasis  is  being  placed  on 
coal  utilization. 

A  significant  amount  of  work  has  previously  been  done  on  the 
direct  utilization  of  coal  in  boilers,  formation  of  soot,  burning  of 
coal  including  temperature  histories,  but  there  is  a  dearth  of 
single  particle  temperature  and  mass  loss  data  reported  in  the 
literature. 

The  goal  of  this  study  was  to  investigate,  record  and  correlate 
the  temperature  profiles  of  individual  graphite  particles  as  they 
burned  and  to  determine  the  mass  lost,  as  a  result  of  laser  heating. 

Preliminary  work  was  done  by  Joseph  Cerv  in  1985  and  is 
documented  in  his  thesis.  To  give  the  reader  some  background  about 
the  importance  of  earlier  work  done  in  this  regard,  I  briefly 
reproduce  Cerv's  introduction  below. 

Theoretically,  the  oxidation  of  carbon  has  been  one  of  the  most 
intensely  considered  topics  in  the  field  of  combustion. 
Experimentally,  the  data  on  carbon  combustion,  especially  dealing 
with  temperature  profiles  and  mass  loss  is  lacking.  The  major 
studies   and   experiments  with   carbon   sphere   combustion  are   as 

follows . 

2 
In  1924,  Nusselt   proposed  his  "shrinking  drop"  model  for  the 

combustion  of  a  solid  sphere  such  as  wood  or  coal.   This  model 
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assumes  that  oxidation  occurs  only  at  the  surface  producing  CO 
and/or  CO  and  that  the  internal  temperature  history  is  governed 
strictly  by  a  one  dimensional  conduction  equation.  When  a  point  on 
the  sphere  reaches  a  characteristics  pyrolysis  temperature,  it  is 
assumed  to  have  burnt.  While  this  model  was  a  start,  it  failed  for 
several  reasons.  It  ignored  temperature  dependence  of  reaction 
rates,  diffusion  limits,  and  ash  effects. 

In  1931,  Burke  and  Schumann  assumed  that  pure  carbon  is 
consumed  due  to  the  reaction  with  CO  at  the  surface  to  form  CO. 
The  CO  thus  formed  reacts  with  the  surrounding  oxygen  at  a  flame 
sheet  to  form  CO  .   The  CO,  thus  formed  at  the  flame  sheet  then 

diffused  back  to  the  surface  to  sustain  the  process. 

3 
In  1951,  Spalding   developed  the  principles  of  liquid  drop 

combustion  which  used  an  evaporative  vapor-phase  diffusion  flame  to 

describe  the  combustion  of  droplets.   While  this  model  adequately 

described  several  solid  and  liquid  combustion  systems,  it  failed  in 

the  case  of  carbon.   Due  to  the  high  sublimation  temperature  of 

carbon  which  is  3800°  K,  it  does  not  have  a  significant  vapor  phase 

concentration  at  realistic  combustion  temperatures. 

In  1958,  Coffin  developed  a  system  for  determining  the  steady 

state  burning  rate  of  solids  utilizing  mass  transfer  coefficients 

rather  than  vapor  phase  diffusion.   Using  this  method,  he  was  able 

to  effectively  model  for  the  cases  in  which  the  mass  transfer 

processes  are  rate-determining. 


The  late  1970 's  brought  expensive  oil  and  renewed  interest  in 
coal  combustion.  This  factor,  coupled  with  an  increase  in  the 
availability  of  computational  facilities,  spawned  several  computer 
model  studies.  Principle  among  the  computer  model  studies  was  the 
one  by  Libby  and  Blake.  They  used  a  computer  model  of  a  single 
coal  particle  that  considered  ash  effects  and  reaction  rates  as  well 
as  diffusion,  established  the  concept  of  extinction  diameter.  The 
extinction  diameter  is  the  particle  size  at  which  a  burning 
particle's  heat  loss  exceeds  its  heat  production.  At  the  extinction 
diameter  the  particle  extinguishes  itself. 

Laurendeau,  in  "Heterogeneous  Kinetics  of  Coal  Char 
Gasification  and  Combustion",  covered  all  aspects  of  the  topic:  coal 
and  char  characteristics,  surface  mechanisms,  particle  reaction 
models,  and  kinetics  of  char  gasification  and  combustion.  Of 
principle  interest  in  this  study  is  the  comparison  of  the  Shrinking 
Drop  Model,  which  strictly  assumes  a  surface  reaction,  to  the 
Progressive  Conversion  Model,  which  considers  pore  effects.  For 
small  particles  on  the  order  of  tens  of  microns,  there  is 
considerable  question  as  to  which  is  the  proper  mechanism. 

Experimentally,  the  literature  is  sparce.  In  the  1930's 
studies  were  done  by  Tu,  Hottel  and  Davis  to  determine  the 
combustion  rate  of  carbon.  Large  spheres  (2.5  cm  in  diameter)  and 
carbon  samples  in  cups  were  heated  in  furnaces,  the  mass  loss  rates 
were  measured  and  the  gas  profiles  around  the  carbon  bodies  were 
sampled.   These  experiments  validated  Burke  and  Schumann's  idea  of 


the  flame  front  for  large  particles.   However,  it  must  be  remembered 
that  this  work  was  conducted  with  large  particles. 

Experimental  work  in  the  area  of  individual  coal  particle 
combustion  slowed  in  the  1940's  as  pulverized  coal  began  being  used 
as  a  fuel  in  boilers.  Interest  shifted  to  furnace  design  with  large 
groups  of  interacting  particles.  Sherman  in  1940  analyzed  the  space 
requirement  for  pulverized  fuel  combustion  using  an  entrained  flow 
of  coal  particles  through  a  furnace.  Even  though  data  was  needed  on 
individual  small  particles,  it  was  not  until  1976  that  work  in  this 

area  was  conducted. 

g 
In  1976,  Ubhayaker  and  Williams   studied  the  burning  rate  of 

small  carbon  particles  (50  urn  -  200  am).   The  particles  were 

suspended  on  a  quartz  fiber  with  grease  and  ignited  with  a  pulse 

from  a  ruby  laser.  As  the  particle  ignited,  it  fell  from  the  quartz 

fiber  and  was  followed  by  a  high  speed  camera.   Additionally,  two 

color  pyrometry  gave  a  temperature  history.   Results  showed  that  for 

40-200  Mm  particles  the  flame  sheet  idea  is  not  valid  and  the  flame 

is   confined   to   the   particle   surface.    Additionally,   particle 

extinction  was  experimentally  verified. 

Since  the  study  by  Ubhayaker,  high  speed  cinematography  and  two 
color  pyrometry  have  been  used  in  several  studies  that  follow 
particles  through  a  furnace.  This  system  gives  burnout  times  and 
temperature  histories  for  50-100  um  particles. 

Real  boilers  use  fuel  sieved  through  200  (74  um)  or  400  (37  um) 
meshes.  This  gives  a  size  distribution  weighted  in  the  20-50  um 
range.   Studies  have  not  been  conducted  on  these  sized  carbon 


particles  because  up  to  now,  there  has  not  been  an  experimental 
method  by  which  they  could  be  measured.  In  the  past  15  years  a  new 
version  of  an  old  device,  the  modified  Millikan  cell,  has  been 
developed  that  can  isolate,  hold,  and  weigh  these  size  particles. 

The  bihyperboidal  Millikan  cell  was  first  developed  in  1970  at 

Aberdeen  Proving  Ground  to  provide  a  method  for  isolating  and 
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observing  20  to  100  urn  aerosol  particles.    In  1981,  Arnold 

developed  the  concept  of  electron  stepping  as  a  method  of  mass  and 

charge  measurement.   This  concept  is  based  upon  the  fact  that  if  a 

charged  particle  lost  an  electron,  the  particle  mass  and  charge  can 

be  determined  from  two  voltage  measurements.   Arnold  applied  this 

concept  to  a  flat  plate  Millikan  cell  and  was  able  to  weigh  picogram 

size  particles.   In  1982,  Philip,  et  al.,  showed  that  the  technique 

of  electron  stepping  could  also  be  used  in  the  bihyperboidal 

Millikan  cell  and  an  extremely  sensitive  mass   isolation  and 

measuring   device  was   available.     It  has   the  capability  of 

capturing  and  holding  a  particle  at  one  position  in  space  and  it  can 

_o        —12 
weigh  particles  in  a  mass  range  from  10   g  to  10    g.   Since  1981, 

the  modified  Millikan  cell  has  been  a  piece  of  equipment  In  search 

of  a  use. 


2.   BACKGROUND 

The   development   of   the   two-color   optical   pyrometer   is 

12 
attributed  to  Hottel  and  Broughton.     They  used  the  pyrometer  to 

measure  true  temperature  and  total  radiation  of  luminous  flames. 

The  two-color  pyrometer  measures  brightness  temperature  at  two 

different,  but  known  wavelengths;  substitution  into  Planck's  law 

yields  two  simultaneous  equations  which  yield  the  true  temperature. 

In  1961,  D.  Rae  used  the  two-color  method  to  determine  the 

temperature  of  frictional  sparks.   Rae  noted  that  sensitivity  of  the 

system  is  greatest  when 


J_  -  1_\    _  2T 


(1) 


and  he  chose  his  photocells  accordingly.  X.  are  the  effective 
wavelengths  of  the  detectors,  and  C  is  a  constant  that  will  be 
defined  later.  Measurements  were  made  in  excess  of  4000°  K.  Rae 
also  used  his  intensity  measurements  to  detect  explosion  of  Titanium 
particles. 

In  1972,  Ayling  and  Smith  burned  pulverized  fuel  particles  in 
a  plug  flow  reactor  that  had  several  ports  for  optical  access.  They 
used  a  single  photomultiplier  with  rotating  filters  to  view  the 
stream  of  particles  passing  by  a  single  port. 

Ubhayaker  and  Williams  were  the  first  to  apply  the  two-color 
pyrometer  to  singly  burning  particles.   They  suspended  electrode 

carbon  particles  in  oil  droplets  on  the  end  of  a  thread  and  ignited 

13 
them  with  a  laser.   House    at  MIT  measured  temperatures  of  coal 


particles  as  they  flowed  past  an  optics  probe  inserted  up  to  the 
bottom  of  a  quartz  laminar  flow  furnace.  Dictor  extended  House's 
work  to  enable  determination  of  complete  temperature  profiles  from 
particle  ignition  to  extinction.  Both  House  and  Dictor  used  a  light 
pipe  to  collect,  split  and  guide  the  light  to  narrow  bandwidth 
detectors.  Signals  from  the  detectors  were  displayed  on  a  memory 
oscilloscope,  photographed  and  measured  by  hand  to  enable 
temperature  calculation.  This  study  is  an  extension  of  the  previous 
work  done  by  Joseph  Cerv  in  the  sense  that  particle  temperature  data 
is  also  incorporated. 


3.   EXPERIMENTAL  DESIGN 

To  enable  temperature  measurements,  a  graphite  particle  must  be 
"isolated";  the  light  from  the  burning  particle  must  be  collected, 
split,  led  to  detectors,  amplified,  measured  and  recorded. 
Additionally,  the  system  requires  calibration,  and  the  raw  data 
needs  to  be  converted  to  actual  temperatures.  The  procedure  to 
accomplish  this  task  is  described  here. 

3.1  Millikan  Cell 

The  equipment  and  technique  for  mass  determination  followed  as 
closely  as  possible  to  that  of  Philip.  The  quadrupole  (Modified 
Millikan  Cell)  shown  in  Figs.  1  and  2  that  suspends  the  particles 
has  surfaces  defined  by  the  following  equations.  The  inner  surfaces 
of  the  endcap  electrodes  are  given  by 


2   r     2 
Z  "  T  +  Zo    •  <2> 


where  Z  »  4mm  and  is  the  vertical  distance  from  the  geometric 
center  of  the  chamber  to  the  endcap  electrodes  and  r  is  the  radial 
distance  from  the  axis  of  symmetry.  The  characteristic  dimension  of 
the  chamber  is  Zq,  half  the  distance  between  the  two  endcap 
electrodes,  measured  along  the  axis  of  symmetry.  The  inner  surface 
of  the  ring  electrodes  is  described  by 


*'£-*.' 
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Fig.  I  Exploded  view  of  Quadrupole 


Fig.  2  Crosa-sectfon    of    Ouadrupola 
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The  quadrupole  with  insulating  spacers  used  in  this  study  was  made 
by  Tempco  Corporation,  Danvers,  Mass. 

The  chamber  consisted  of  a  central  ring  electrode  of  aluminum 
and  two  teflon  rings  on  either  side.  The  two  aluminum  endcap 
electrodes  were  mounted  in  the  teflon  rings.  Ten  (1/8")  access 
ports  were  drilled  into  the  teflon  ring  and  endcap  electrodes.  All 
access  ports  could  be  sealed  by  flat  glass  windows  when  not  in  use. 
All  outside  metal  surfaces  on  the  quadrupole  were  covered  with 
teflon  shields  to  prevent  electrical  shock. 

The  circuitry  to  produce  the  AC  voltage  (with  DC  bias)  on  the 

ring  electrode  and  the  DC  voltage  across  the  endcaps  is  given  in 

Fig.   3.   The  capacitor-resistor  network  suppresses   induced  AC 

voltages  on  the  endcap  electrode  resulting  from  the  AC  voltage  on 

the  ring  electrode.   The  regulated  dc  power  supply  to  produce  V 

across  the  end  cap  consists  of  a  Pacific  Instruments  7104  PC  card. 

This  card  is  capable  of  producing  0-2000  V  across  the  end  caps  with 

a  resolution  of  0.0005  V.   The  power  transformers  were  Essex  Stancor 

P-8151  2400  VAC  CRP  transformers.   The  power  transformers  were 

controlled  by  a  Powerstat  T/16  nongrounded  variable  transformer. 

The  ring  electrode  system  provided  an  rms  ac  voltage  of  0  to  4800  V 

at  60  Hz  with  a  resolution  of  about  50  V.   Measurements  of  V  were 

g 

made  with  a  Keithly  192  digital  multimeter  capable  of  reading  up  to 
6}  digits.  The  meter  was  protected  by  a  30  A  fuse  on  the  negative 
lead. 
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I  10  V    60  Hz. 


Pig.  3     Circuitry    for    Quadrupole 
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3.2  MASS  MEASUREMENT  AND  ELECTRON  STEPPING 

The  schematic  of  Che  apparatus  for  mass  measurement  and 
combustion  optics  is  shown  in  Figs.  4  and  5,  respectively. 

A  charged  particle  is  levitated  between  a  set  of  hyperbolic 
electrodes,  the  schematic  diagram  of  which  is  shown  in  Fig.  2. 
There  are  two  endcap  electrodes  between  which  a  DC  potential  is 
applied,  as  well  as  a  ring  electrode  to  which  an  AC  voltage  is 
applied.  The  DC  potential  across  the  endcaps  balances  the  particle 
against  gravity,  while  the  AC  voltage  on  the  ring  provides  for  the 
focussing  of  the  particle  towards  the  center  of  the  chamber.  The 
particle  is  illuminated  by  a  HeNe  laser  beam  that  traverses  the 
chamber  through  ports  drilled  into  the  center  of  the  endcap 
electrodes.  As  shown,  the  scattered  light  from  the  particle  is 
observed  through  ports  in  the  endcap  electrode  at  23°  to  the  HeNe 
beam. 

At  the  geometric  center  of  the  device  there  is  no  AC  field. 
Therefore  when  the  electric  force  produced  by  the  DC  field  just 
balances  the  gravitational  force  on  the  particle,  the  center  of  the 
chamber  becomes  a  stable  equilibrium  point  for  the  particle.  Stated 
differently,  whenever  a  balanced  particle  drifts  away  from  the 
center  it  "feels"  a  net  force  towards  the  center  of  the  chamber. 

Experimentally,  the  mass  is  measured  as  follows.   The  ring 

voltage  is  initially  adjusted  to  about  500  V.   The  voltage  across 

the  endcaps  was  set  to  its  minimum  value,  typically  about  0.015  V. 

Particles  were  injected  into  the  chamber  with  a  hypodermic  needle. 

Due  to  the  applied  AC  field  in  the  ring  electrode,  this  often  leaves 

several  charged  particles  stably  oscillating  in  the  quadrupole. 
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Then  a  particle  that  responded  suitably  to  the  balancing  DC  voltage 
was  balanced  at  the  center  of  the  chamber.  A  stationary  particle 
was  assumed  to  be  balanced  since  any  departure  from  the  center  of 
the  chamber  resulted  in  oscillations  of  the  particle.  The  ring  AC 
voltage  was  then  lowered  for  a  few  seconds.  This  usually  allowed 
the  unbalanced  particles  to  fall  out  of  the  chamber  while  the 
balanced  particle  remained  near  the  chamber  center.  The  AC  voltage 
was  then  increased  again.  The  remaining  particle  was  more  carefully 
balanced  by  adjusting  V  and  then  lowering  the  AC  voltage  for  a 
couple  of  seconds.  If  the  particle  drifted  from  the  chamber  center, 
it  was  clearly  unbalanced.  The  ring  voltage  was  then  increased 
again  and  V  further  adjusted,  after  which  the  procedure  was 
repeated,  until,  on  lowering  and  raising  the  AC  voltage  the  particle 
remained  stable,  which  is  recorded.  The  UV  light  was  then  turned  on 
for  an  instant.  Electron  loss  was  verified  by  lowering  the  AC 
voltage  and  observing  particle  drift.  The  particle  was  then 
rebalanced  at  the  chamber  center  by  the  technique  already  described 
and  Vg  recorded  again.  This  procedure  is  repeated  a  number  of  times 
and  from  several  measurements  of  V  the  mass  and  charge  of  the 
particle  can  be  calculated. 

3.2.1  Mass  Calculation 

For  a  stationary  particle  at  the  center  of  the  quadrupole,  a 
simple  force  balance  describes  the  mass  of  the  particle  as 

mg  -  qE 
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where  E  is  the  electric  field.   But  for  the  modified  Millikan  cell 

its  value  is  C  V  /Z  . 
g  o 

q  is  particle  charge 

C  is  the  geometric  constant  for  the  chamber  equal  to  0.4 

g  is  the  gravitational  constant 

Z  is  the  characteristic  dimension  of  the  quadrupole  Z   ■  4mm 
and  V  is  the  potential  difference  across  the  endcaps . 

The   charge  on   the  particle   equals   the   number   of   excess 
electrons  times  the  charge  on  an  electron,  so 

Cne  V 


m  =  -^   .  (4) 


Collecting  constants 


where     K  =  Ce/gZ  =  1.6356  x  10"15  g/V. 


If  vn+i  ami  v  +i   represent  the  balancing  voltages  for  a 
particle  containing  n  and  n+1  excess  electrons,  respectively,  then  by 
solving  the  previous  equation  at  each  balancing  voltage  the  particle 
mass  is  given  by 


K  Vn  Vn+1 

V  I~7     '  ^ 

n    n+1 


For  relatively  large  particles  (10~  )  with  large  numbers  of 
excess  electrons  (10  )  over  a  small  range  of  voltage  (lv)  this  can 
be  approximated  as 
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(6) 


AV  for  a  single  electron 

The  problem  is  now  reduced  to  finding  the  voltage  change  for  a 
single  electron. 

3.2.2  Measurement  of  V 
g 

There  are  two  methods  used  to  measure  V  .   The  first  method  is 

g 

to  capture  a  particle  in  the  center  of  the  chamber  and  manually 
adjust  V  until  the  particle  does  not  oscillate.  The  AC  voltage  is 
turned  down  and  the  direction  of  the  drift  is  noted.  The  AC  field 
is  turned  alternately  up  and  down  and  V  is  adjusted  until  the 
particle  drifts  only  horizontally.  The  accuracy  of  this  method  is 
limited  by  the  ability  of  the  experimentalist  to  resolve  the  change 
in  the  direction  of  drift  with  a  change  in  a  balancing  voltage. 

The  second  method  is  to  record  the  voltage  from  the  position  of 
the  particle  as  the  charges  are  removed.  The  position  of  the 
particle  is  fixed  with  respect  to  the  cross  wires.  This  method  was 
found  to  be  less  accurate.  The  particle  position  is  adjusted  by 
manipulating  the  AC  and  DC  fields  in  such  a  manner  that  it  "touches" 
the  x  and  y  cross  wires.  The  balancing  voltage  is  recorded  (V.). 
After  irradiating  with  UV  light  the  particle  drifts  to  a  new 
position.  By  adjusting  the  DC  voltage  (V  )  the  particle  is  brought 
back  to  it's  former  position.  A  series  of  measurements  are  made  to 
evaluate  AV. 
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3.2.3  Particle  Sizing 

Illumination  was  provided  by  Spectra  Physics  Model  102-3  HeNe 
laser.  The  laser  was  mounted  on  a  single  axis  translational  stage 
and  attached  to  a  lead  brick.  This  mounting  coupled  with  a  beam 
steerer  allowed  the  laser  beam  to  be  directed  into  the  chamber  along 
the  asymptote  between  the  electrodes.  The  particles  were  observed 
and  sized  through  a  flat  glass  window  at  23°  to  the  illuminating 
beam  using  a  Beck  telescope  with  an  image  distance  of  200  mm,  a  50 
mm  objective,  and  a  10X  Beck  micrometer  eyepiece.  Sizing  was 
accomplished  by  sweeping  the  movable  reticle  across  the  particle 
image.  Adjusting  the  AC  field  to  align  one  edge  of  the  particle 
with  the  stationary  stadia  line  gives  the  best  results.  This 
optical  system  allowed  a  resolution  of  2.5  urn  in  sizing  particles. 
The  UV  light  for  electron  stepping  was  provided  by  an  Oriel  6035  low 
pressure  mercury  calibration  lamp.  The  lamp  was  mounted  in  one  of 
the  ten  chamber  access  ports. 

3.3   RADIATION  THEORY  AND  TWO-COLOR  PYROMETRY 

The  task  of  obtaining  time  temperature  profiles  of  singly 
burning  graphite  particle  is  not  a  simple  one.  Temperature  might  be 
low,  the  total  duration  of  the  burnout  may  be  a  few  milliseconds, 
and  particle  sizes  are  on  the  order  of  microns.  The  nature  of  this 
project  precluded  the  use  of  any  nonoptical  means  for  measuring 
surface  temperature. 

The  light  intensity  of  a  burning  graphite  particle  at  any 
wavelength  is  proportional  to  the  particle's  surface  area  and  a 
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strong   function   of   temperature   and   emissivity.    The   radiant 
intensity  is  given  by  Planck's  law  (Fig.  6). 

2Trhc      e. 


[exp  (iff. 


h  =  — 5 —  lexp  ISri  "  l  C7) 


where  c  is  speed  of  light,  h  is  Planck's  constant,  K  is  Boltzmann's 
constant,  T  is  particle  temperature,  A  is  the  wavelength  and  e.  is 

the  emissivity  at  that  wavelength.   The  emissivity  is  a  function  of 

2 
temperature.   2nhc   is  Planck's  first  radiation  constant,  C  ,  and 

-12      -2 
equals  3.74  x  10    W  cm   .   Similarly,  hc/K  is  Planck's  second 

radiation  constant,  C  ,  and  equals  1.44  cm  °K.   For  exp  C./XT  >>  1, 

i.e.   for  a  maximum  wavelength  of  observation  of  900  nm  and 

temperature  below  3500°  K,  we  can  use  Wein's  approximation  and 

disregard  the  unity  term  on  the  right  hand  side  of  Eq.  (7).   The 

radiative  flux  at  X  is  now  given  by 


-  exp  [wj 


-C„ 

(8) 


Particle  temperature  cannot  be  directly  determined  using  Eq.  (8). 

The  emissivity  is  unknown.   In  two-color  pyrometry  one  is  concerned 

not  with  absolute  intensity  at  one  wavelength  but  rather  with  the 

ratio  of  two  intensities  from  two  different  wavelengths  A.  and  X  . 

s     A      B 

Assuming  that  emissivity  is  not  a  function  of  wavelength  (grey 
body)  ,  the  ratio  e,/eB  becomes  unity  and  does  not  appear  in  the 
ratio  of  the  emitted  power  at  two  different  wavelengths, 


03       1  5        10  50       W* 

Wavelength  (urn) 

Fig.  fo    The  Planck  distribution  law;  spectral  radiance  of  blackbody  radiation  as  a  function  of 
temperature  and  wavelength. 
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(A+Xl   EXA  « 
R(T,X)  =  -r .  (9) 

LB+Al  E» dA 

>AB-X, 

The  measured  intensity  at  each  wavelength  will  depend  on  system 
parameters  such  as  the  amount  of  light  carried  to  the  detectors  and 
the  detector  efficiency.  The  latter  is  often  dependent  on 
wavelength.  If  the  bandwidths  are  narrow,  we  can  assume  the  average 
efficiency  of  the  detectors  within  a  bandwidth  to  be  that  of  the 
peak  wavelength.   Equation  (9)  can  be  approximated  as 

-  _  h       \,A  "p.A  XB  AXA     fC2  fl    111 

i 5 exp  f~  u  "  T   •   (10) 

where  I.  are  the  measured  intensities,  AX.  are  the  bandwidths  used 
in  detection,  nL  ,  are  the  efficiencies  for  carrying  light  to  the 
detectors  and  rtD  ,  are  the  detector  efficiencies  at  each  wavelength. 
All  of  the  coefficients  for  the  exponential  term  in  the  above 
equation  can  be  combined  into  one  calibration  constant  leaving 

=.-.p  P  !-- 


P  P  -  *-l 


Equation  (11)  is  graphed  in  Fig.  7.   Inverting  Eq.  (11)  one  finds 

n        fl      1 


XA        h 


In  R  -  In  C   - 

This  approximation  is  valid  for  narrow  bandwidths  only. 
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3.3.1   Effect  of  Bandwidth  on  Ratio  of  Intensities 

The  bandwidth  of  the  filters  used  was  approximately  100  nm. 
This  enabled  us  to  obtain  a  larger  signal  and  hence  measure  lower 
temperatures.  Detailed  calculations  (Sec.  3.3.2)  were  made  to  see 
the  effect  of  different  bandwidths  on  the  ratio  of  the  intensities 
at  two  chosen  wavelengths.  Rather  than  using  a  point  to  point 
ratio,  integrated  "bands"  beneath  the  radiation  curve  must  be  used. 
Corrections  were  made  for  filter  transmission  (part  of  n,  )  and 
detector  efficiency  r\   )    so  Eq.  (9)  is 

XA+AAA/2 

R  =  c  **& " — 7^     (13) 

\,BU)  VB(X)  XB5  eXp  P^Jj  AXB 


2 


w 


The  functionalities  of  n.  .  and  n« .  have  to  be  determined.   If  the 
Li      Di 

bandwidths  become  too  wide  and  being  to  overlap  one  another,  the 
sensitivity  of  IA/IR  to  temperature  diminishes.  The  calibration 
constant  (C)  now  arises  from  the  geometry  of  the  system. 


3.3.2  Effect  of  Bandwidth  on  Temperature  Measurement 

From  Eq.  (11)  it  is  evident  that  the  ratio  of  the  two 
intensities  I  and  I  gives  the  temperature  of  the  particle.  Using 
Planck's  law  the  ratio  of  intensities  at  the  two  wavelengths, 
namely,  900  and  750  nm  were  calculated  for  different  temperatures  as 
a  function  of  the  filter  bandwidths.  The  results  are  listed  in 
Table  1  and  graphed  in  Fig.  8. 
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TABLE  1 


WAVELENGTHS     BANDWIDTH    TEMPERATURE     RATIO  I./I„ 

A  B 


in  nm         in  run        in  K 


900     750        10         2000  1.851 

2500  1.384 


3000  1.142 


900     750        40         2000  1.850 

2500  1.385 


3000  1.143 


900     750        80         2000  1.856 

2500  1.385 


3000  1.144 


900     750        120         2000  1.838 

2500  1.384 


3000  1.145 


900     750        200         2000  1.815 

2500  1.383 


3000  1.149 


900     750       400         2000  1.717 

2500  1.372 


3000  1.165 
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FIGURE  9a 

Intensity  trace  for  a  lignt  bulb  using  a  10  nm  filter.  Tne 
experimental  set-up  consisted  of  a  "cnopper"  in  between  an  AC 
lignt  bulb  and  tne  pyrometer.   Tne  "pulse"  seen  in  tne  figure 
is  approximately  300  Us  wide  in  tne  duration  for  wnicn  tne 
pyrometer  could  "see"  tne  lignt  bulb. 


FIGURE  9b 

Intensity  trace  for  a  lignt  bulb  using  a  100  nm  filter.  Time  scale 
500  us/division. 
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The  analysis  shows  that  the  ratio  of  the  intensities  at  the  two 
wavelengths  does  not  change  significantly  even  if  the  bandwidth  of 
the  filters  is  varied  from  10  to  100  nm,  indicating  that  the 
temperature  as  measured  from  the  ratios  would  be  the  same,  but  on 
the  other  hand  the  signal  detected  with  a  100  nm  filter  would  be  10 
times  that  of  a  10  nm  bandwidth  filter  (Figs.  9a  and  9b).  Hence  the 
signal  to  noise  is  greatly  improved  by  using  a  100  nm  bandwidth 
filter  without  sacrificing  the  ratio,  thereby  giving  a  better 
accuracy  in  measuring  the  temperature. 

3.4  ACCESS  TO  BURNING  PARTICLE 

A  three  foot  bifurcated  fiber  optic  cable  manufactured  by 
Ealing  was  used  to  transmit  the  light  signal  from  the  particle  to 
the  pyrometer.  A  schematic  of  this  set  up  is  shown  in  Fig.  10.  The 
advantages  of  using  the  fiber  optic  cable  are  manyfold;  namely: 

i)  Radiation  from  the  burning  particle  is  received  only  a  short 
distance  from  the  particle  thereby  eliminating  the  need  for 
optical  focussing  inside  the  Millikan  Cell. 
Ii)  Reduction  in  noise  level  of  the  pyrometer  since  the  AC 
interference  of  the  ring  electrode  is  minimized. 
lii)  Permitted  the  design  of  the  filter  detector  assembly  and  the 
signal  receiving  end  of  the  system  as  separate  units.  The 
resulting  versatility  of  the  pyrometer  has  proved  to  be  very 
advantageous . 
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3.4.1   Particle  Simulation 

Here  I  would  demonstrate  the  ability  of  our  pyrometer  to 
measure  temperature  of  a  "typical"  particle.  Consider  a  graphite 
particle  of  radius  (r) ,  burning  at  a  temperature  T.  The  total  black 
body  power  emitted  by  the  particle  at  wavelength  X  in  the  interval 
dX   is  given  by  the  expression 


27rhc^   dX 

0  '^] — n^~r-T  • 


X    [eXpllKT|-! 

This  power  is  radiated  in  all  directions  by  the  particle,  a 
fraction  of  which  is  incident  upon  the  fiber  optics.  The  power 
incident  on  the  fiber  is 

P  -  T.P  x  A  x  S 
o 

where     T.P  is  total  power  emitted  by  the  particle 
A  is  the  surface  area  of  the  particle 

S  is  the  solid  angle  the  fiber  subtends  on  the  particle 
Considering   losses   associated   with   the   fiber,   filters, 
photodiodes  and  so  forth  (Fig.  11)  the  voltage  signal  detected  by 
the  oscilloscope  is 


V  =  T.P  x  A  x  S   x  L.  x  L.  x  L,  x  S  x  R. 
o     1    2    3    e    f 


where      P  is  power  incident  on  fiber 

L  is  the  loss  due  to  bifurcated  fiber  (50%) 
L9  is  the  transmission  loss  in  fiber  (507.) 
L  is  the  transmission  loss  in  filter  (50%) 
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Se  is  the  sensitivity  of  the  photodiode  (.3  A/W) 
R,  is  the  feedback  resistor  of  the  OPAMP  (200  MS}) 
For  a  typical  "particle"  of  radius  20  urn,  burning  at  a 
temperature  of  1300°  K,  the  total  power  emitted  at  900  nm  (100  nm 
interval)  is  298.4  W/m2. 

Surface  area  of  the  particle       5.026  x  10   m2 

Solid  angle  subtended  by  fiber     5.625  x  10 
situated  1  cm  from  the  particle 

Loss  associated  with  bifurcation   0.50 
of  fiber 

Loss  associated  with  fiber  length  0.50 

Transmission  loss  with  filters     0.50 

Sensitivity  of  photodiode  at       0.3  A/W 
900  nm 

Feedback  resistance  of  OPAMP       200  Mfi 

Substituting  the  above  values  in  the  previous  equation  gives  a 
voltage  signal  of  60  mv.  The  noise  level  of  the  amplifier  was 
10  mv. 

3.5   PYROMETER  DESIGN 
3.5.1  Choice  of  Detectors 

The  temperature  of  the  burning  graphite  particles  was  expected 
to  be  in  the  range  of  1300  -  3000°K.  A  plot  of  the  emitted  energy 
vs  wavelength  for  a  black  body  source  is  given  in  Fig.  6.  From  the 
plot,  in  the  temperature  range  of  interest,  the  amount  of  radiant 
energy  emitted  at  wavelengths  below  700  nm  is  very  small.   In  view 


33 


of  this  it  was  decided  to  use  photodiodes  rather  than  PMT '  s  which 

have  a  better  response  in  the  infrared. 

A  second  consideration  in  the  choice  of  wavebands  arises  from 

the  grey  body  assumption  inherent  in  the  determination  of  the 

temperature  from  the  signal  ratios  for  unknown  e,   and  e,  .   The 

assumption  that  e^  =  eAB  is  best  for  AA  close  to  A     A  third 

consideration  is  that,  for  a  good  precision  in  determining  T  from  R, 

dR/dT  should  be  large.   dR/dT  is  greatest  in  the  visible  and  near 

infrared  (A  close  to  A  ). 
A  B 

The  chosen  bands  A   =  900  nm,  A„  =  750  nm  and  AA,  =  AA   = 
A  B  A      B 

100  +  20  nm  met  all  the  requirements.  Bandpass  interference  filters 
matching  these  specifications  were  obtained  from  Andover 
Corporation.  The  peak  transmission  of  these  filters  is  about 
65-70Z.  The  photodiodes  chosen  were  S874-5K  silicon  photocells 
manufactured  by  Hamamtsu  Corporation.  The  sensitivity  is 
approximately  0.3  A/W  at  these  two  wavelengths  and  the  active 
photosensitive  area  is  5.6  mm2. 

3.5.2  Signal  Processing 

Two  amplifiers  with  a  gain  of  2  x  10  were  built  to  amplify  the 

detected  signals.   RCA  SK9378  Bimos  OPAMP  were  chosen  for  their  fast 

response  and  low  noise.   The  schematic  of  the  current  amplifier  is 

shown  in  Fig.  12.   Because  of  the  small  currents  being  measured 

-9 
(10    Amp),  shielding  considerations  were  important  for  the  input 

signals  and  amplifiers. 
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3.5.3  Shielding  Consideration 

Even  though  the  amplifier  circuit  looks  quite  simple 
(consisting  of  an  OPAMP,  resistor  and  capacitor),  shielding 
considerations  are  vital  in  designing  circuits  dealing  with  low 
level  currents. 

Once  the  OPAMP  (SK  9378)  was  chosen,  the  next  step  was  to 
optimize  circuit  performance.  Because  the  amplifiers  extraordinary 
sensitivity  makes  high  performance  applications,  error 
considerations  which  are  insignificant  in  less  sophisticated 
applications  now  become  the  limiting  factor  in  circuit  operation. 

Hard  wiring  with  good  quality  insulation  (teflon)  was  preferred 
to  printed  circuit  wiring.  The  amplifier  was  plugged  into  a  socket 
with  teflon  insulators  rather  than  soldered  directly  onto  a  printed 
circuit  board.  Insulated  wires  connected  input  terminals  to 
external  components.  Further,  all  wiring  was  short  and  supported  in 
such  a  way  that  movement  due  to  vibration  is  limited  as  much  as 
possible,  because  motion  of  input  wiring  will  cause  corresponding 
variations  In  parasitic  capacitance. 

Coaxial  cable  was  used  wherever  possible  to  minimize  RF  pickup. 
Metal  film  resistors  were  chosen  for  their  low  noise,  good 
temperature  stability,  low  shunt  capacitance  and  relatively  low 
cost.  Polystyrene  capacitors  were  used  for  decreased  leakage 
resistance.  Finally,  the  amplifier  and  photodiodes  were  enclosed  in 
an  aluminum  container  to  shield  the  60  Hz. 


36 


3.6  Calibration  of  Pyrometer 

Calculation  of  temperature  information  from  the  pyrometer 
outputs  required  determining  the  constant  in  Eq.  (111.  Calibration 
was  performed  using  a  Sylvania  DXW  1000  W  bulb  operating  at  a 
standard  temperature  of  3200  +  50 °K.  Information  about  the  quality 
of  this  lamp  and  its  use  as  a  source  of  spectral  irradiation  is 
provided  by  the  Infrared  Handbook.  7  For  actual  calibration  runs 
with  the  DXW  lamp,  the  pyrometer  was  shielded  from  stray  light  rays 
and  the  two  base  lines  (900  &  750  nm)  were  made  to  coincide  at  a  set 
line  on  the  dual  trace  scope.  The  DXW  lamp  was  placed  at  a  fixed 
distance  from  the  pyrometer  and  the  response  photographed  on  a 
polaroid  film.  To  prevent  reflected  light  from  entering  the 
pyrometer,  the  optical  fiber  was  encased  in  a  black  casing.  This 
proved  effective  since  the  ratio  of  intensities  did  not  vary 
significantly  even  if  the  lamp  was  moved  from  the  pyrometer.  The 
calibration  constant  can  be  calculated  from  Eq.  (12). 

f. 


1 J_ 

X„  "  X. 


Jin  R  -  Jin  C 


Substituting  the  value  of  the  known  quantities  Viz.  T  ■  3200°K, 
C£  •  1.44  cm°K,  XA  =  900  nm,  Xg  =  750  nm  and  the  ratio  of  the 
intensities  R  calculated  from  the  polaroid  photograph,  the  constant 
C  can  be  evaluated.  The  constant  C  varied  by  8%  In  10  days  as  a 
function  of  time  and  hence  had  to  be  determined  before  each  run. 
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3.7   PYROMETER  RESPONSE  TIME 

Figure  13  describes  the  schematic  for  determining  the  pyrometer 
response  time.  The  mechanical  chopper  consisted  of  a  circular  disk 
driven  by  a  motor.  A  small  aperture  (at  a  distance  of  5.5  cm  from 
the  center  of  the  disk)  was  made  to  allow  "chopping"  of  the  HeNe 
laser  beam.  The  laser  beam  was  focussed  (convex  lens)  on  the  plane 
of  the  disk.  The  time  period  of  revolution  of  the  chopper  was 
determined  from  the  oscilloscope  trace.  One  can  find  the  linear 
velocity  at  the  "aperture"  using  the  simple  relation, 


2-f   ■  (15) 


where     V  is  the  linear  velocity  at  the  "aperture". 

R  is  the  distance  between  the  "aperture"  and 
center  of  the  disk.   (5.5  cm). 

T  is  the  time  period  of  revolution  (24  ms) . 


Substituting  the  above  values  in  Eq .   (15)  we  get  V  as 

3 
1.44  x  10  cm/sec.   The  beam  radius  of  the  laser  at  the  focal  point 

can  be  obtained  from  the  relation 


"i-rr-  (16) 

o 

where     X  is  the  wavelength  of  the  laser  beam  (0.63  urn), 
f  is  the  focal  length  of  the  lens  (5.7  cm), 
u  is  the  beam  radius  of  the  unfocussed  beam  (1  mm). 


Hence  the  diameter  of  the  focussed  beam  (2ai  )  is  2.3  x  10    cm. 
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FIGURE  Ha 

Pyrometer  output  for  response  time. 

Time  scale  100  us/div. 

Offset  voltage  of  OPAMP  20  mV. 


FIGURE  14b 

Pyrometer  output  for  response  time. 

Time  scale  100  us/div. 

Offset  voltage  of  OPAMP  800  mV. 
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This  gives  a  theoretical  chopping  time  (2id . /V)  of  15  us.   The 
response  time  of  the  pyrometer  as  shown  in  Fig.  (14b)  is  200  us. 

3.8   EXPERIMENTAL  PROCEDURE 

This  section  details  the  sequence  of  operation  for  data 
collection.  At  the  beginning  of  each  run,  the  calibration  constant 
C  was  determined  as  described  in  Section  3.6.  The  following 
procedure  was  followed  for  each  set  of  experimental  conditions. 

i)   A  single  graphite  particle  "trapped",  sized  and  massed  as 
described  in  Section  3.2. 
ii)   By  manipulating  the  DC  field,  the  particle  was  made  to  pass 
through  the  Ar  ion  laser  beam  until  a  trace  was  obtained  on 
the  memory  oscilloscope, 
iii)   The  particle  was  massed  and  the  temperature  was  measured  from 
the  ratio  of  intensities. 
iv)   Repeat  conditions  starting  at  (i)  for  a  "fresh"  particle  or 
same  particle  if  still  present  in  the  field  of  view. 


4.   RESULTS  AND  DISCUSSION 

A  variety  of  information  can  be  derived  from  the  pyrometer 
outputs.  The  ratio  of  the  outputs  at  two  different  wavelengths 
gives  the  temperature.  A  single  pyrometer  output  depends  upon  the 
particle  size  and  the  ambient  environment  around  it. 

4.1  PYROMETER  OUTPUT  FOR  A  LIGHT  BULB 

A  standard  75  W  light  bulb  operating  at  110  v  and  60  Hz  was 
used  to  test  the  response  of  the  pyrometer.  The  experimental  set-up 
is  similar  to  Fig.  10  except  that  the  HeNe  laser  was  replaced  by  the 
light  bulb.  The  pyrometer  output  and  the  temperature  profile  are 
given  in  Fig.  (15).  The  reader  is  encouraged  to  find  the  AC 
frequency  from  the  pyrometer  output. 

4.2  SINGLE  PARTICLE  TRACES 

All  of  the  raw  data  in  this  study  was  in  the  form  of  two  traces 
on  the  screen  of  a  memory  oscilloscope  (900  nm  at  the  top  and  750  nm 
at  the  bottom).  Unless  otherwise  stated  the  time  scale  is  500us/div 
and  voltage  scale  is  lOOmv/div  on  all  the  traces.  Polaroid  pictures 
of  the  traces  used  to  find  time-temperature  histories  are  listed  in 
the  appendix. 

Figure  (16)  is  a  typical  polaroid  trace.  The  signals  are 
inverted  and  appear  here  as  they  did  on  the  oscilloscope.  Time 
increases  from  left  to  right.  The  time-temperature  analysis  was 
computed  as  given  in  Table  2. 


43 


FIGURE  15a 

Intensity  trace  of  a  50W  AC  lignt  bulb.   Tne  top  and  bottom  traces 
are  using  900  and  750  nm  filters.   Time  scale  5  ms/div. 


FIGURE  15b 
Temperature  profile  for  Fig.  15a. 
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TABLE  2 


0.95/0.5 

1593 

1.9/1.1 

1679 

3.4/2.1 

1743 

3.95/2.15 

1622 

3.6/1.8 

1550 

2.7/1.3 

1519 

1.9/0.7 

1335 

1.4/0.5 

1317 

1.2/0.37 

1235 

Time  in  Ratio  of  Intensity  at       Temperature 

Microseconds        900  and  750  nm  (cm/cm)       in  Kelvin 


500 
700 
900 
1100 
1300 
1500 
1700 
1800 
1900 


The  temperatures  were  calculated  from  the  ratios  using  Eq.  12. 
The  value  of  the  calibration  constant  C  was  found  by  substituting 
the  ratio  of  intensities  at  900  and  750  nm  obtained  from  the 
pyrometer  response  to  the  1000W  DXW  bulb  operating  at  3200  ±  50°K. 
The  ratio  in  this  case  was  0.75  thereby  giving  a  value  of  ln(C)  ■ 
1.208.  Temperature  data  for  points  from  0  to  500  us  were  not 
computed  due  to  very  low  signals  at  750nm  (few  millivolts).  An 
error  of  0.1cm  in  determining  the  height  (I  )  in  this  region  could 
result  in  an  error  of  200°K  in  the  temperature  (Fig.  7). 

Oscillating  intensities  can  be  seen  in  Fig.  17.   There  exist  at 
least  four  possible  explanations  for  the  oscillations: 
i)   temperature  is  fluctuating 
ii)   the  particle  rotates  exposing  different  surface  areas  to  the 
detector 
iii)   the  particle  is  moving  in  and  out  of  the  field  of  view 
iv)   "ringing"  in  the  amplifiers. 
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FIGURE  16a 

Pyrometer  output  for  a  30um  grapnite 
particle  neated  by  a  laser. 

Laser  Power  0.32  W 

Ambient  environment  OXYGEN 
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FIG.   17 
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FIG.   18 


52 


FIG.   19A 
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The  profiles  from  several  traces  were  analyzed  to  see  If  (i) 
was  true.  No  exact  correlation  was  seen  between  absolute  intensity 
and  temperature  even  though  temperature  did  fluctuate.  For  traces 
showing  oscillations  and  a  definite  rise/fall  behavior  in  absolute 
intensity,  a  curve  was  drawn  through  the  middle  of  each  trace. 
Measurements  were  made  using  this  curve  (Fig.  17) .  The  error 
introduced  by  this  method  was  quite  large.  Figure  18  is  included  to 
give  an  example  of  a  "suspicious"  trace.  It  is  quite  possible  that 
a  second  particle  (obtained  due  to  fragmentation  of  the  "parent" 
particle)  ignited  before  the  first  had  extinguished,  and  the  traces 
from  both  appear  as  one  on  the  oscilloscope.  This  explanation  was 
confirmed  in  some  cases  by  visual  observations.  Possibility  (iii) 
was  ruled  out  since  the  diameter  of  the  optical  fiber  was  3mm 
OOOOum)  and  the  particle  oscillations  were  confined  to  800um  only. 
The  author  investigated  (iv)  further  since  it  was  the  prime 
"suspect".  The  pyrometer  response-time  experiment  was  performed, 
with  the  time  scale  at  lOOus/div.  The  inquisitive  reader  might 
notice  (Fig.  Ha)  "fine"  oscillations  occurring  at  an  interval  of 
about  300  us  or  3  KHz.  This  was  the  "ringing"  frequency  found  in 
most  of  the  traces  showing  oscillatory  behavior.  Efforts  to 
smoothen  these  oscillations  by  changing  the  offset  voltage  of  the 
OPAMP's  to  several  hundred  millivolts  were  futile. 

In  most  cases  10  or  more  data  points  were  determined  for  each 
trace.  Sufficient  data  was  collected  to  justify  the  smooth  curve 
approximation,  but  the  main  reason  for  providing  the 
time-temperature  graphs  is  to  illustrate  how  temperature  can  vary 
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within  a  single  trace.  For  a  majority  of  the  traces  there  appears 
to  be  a  rise-fall  behavior  in  temperature  as  well  as  absolute 
intensity  though  there  are  traces  where  temperature  rises  as  the 
absolute  intensity  falls.  Figures  19a  and  19b  illustrate  this  point 
and  this  might  well  be  an  experimental  "artifact".  In  all  the 
time-temperature  graphs,  the  data  points  are  connected  by  a  straight 
line.  This  could  be  somewhat  inappropriate  and  misleading,  but  is 
the  clearest  and  best  representation  of  the  data  available. 

4.3  MEAN  PARTICLE  BEHAVIOR 

Several  runs  were  made  to  test  the  effects  of  ambient 
environment  (oxygen,  nitrogen  and  air)  and  particle  size.  A 
summary  of  the  data  is  presented  in  Table  3  and  Figures  20  through 
28. 

TABLE  3 
Ambient  environment  AIR 
Laser  power         0.30W 
Laser  beam  diameter   22  urn 


Particle  Size     Peak  Temperature     Particle  Heating 
in  Microns         in  Kelvin        time  in  us  (FWHM) 


20  2137  1100 

22  1772  1400 

22  2187  900 

27  1752  3500 

27  2573  800 

33  1868  600 
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Ambient  environment    NITROGEN 
Laser  power  0.30-0.32W 

Laser  beam  diameter    22  urn 


Particle  Size 
in  Microns 


Peak  Temperature 
in  Kelvin 


Particle  Heating 
time  in  us  (FWHM1 


25 

1738 

25 

1740 

Above  particle 

heated  again 

1679 

27 

1847 

27 

1758 

30 

1427 

Above  particle 

heated  again 

1629 

30 

1654 

Above  particle 
heated  again 


1834 


1300 
1000 

2100 

500 

1200 

900 

900 
1200 

1200 


Ambient  environment     OXYGEN 
Laser  power  0.34W/0.322* 

Laser  beam  diameter     22  um. 


Particle  Size 
in  Microns 

Peak  Temperature 
in  Kelvin 

Particle  Heating 
time  in  us  (FWHM1 

15 

2059 

1200 

27 

1883 

2100 

27 

2648 

400 
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Particle  Size     Peak  Temperature 
in  Microns         in  Kelvin 


Particle  Heating 
time  in  us  (FWHM) 


30 

1624 

30 

1762 

20* 

2335 

20* 

2115 

29* 

2035 

Above  particle 

heated  again 

1790 

30* 

1641 

Heated  again 

1743 

Heated  again 

1736 

32* 

1538 

34* 

1724 

Heated  again 

1842 

>2800 

1300 

300 

900 

1300 

>2200 

1300 

1100 

900 

700 

2000 

>3500 


Neglecting  any  of  the  smaller  errors  that  may  have  arisen  in 
the  course  of  collecting  and  analyzing  data,  I  would  estimate  an 
overall  error  in  temperature  of  50  to  150°K. 

Despite  the  unusual  behavior  of  the  individual  graphite,  some 
consistent  trends  in  average  temperature  profiles  have  been  seen 
(Table  3  and  Fig.  28).  The  most  notable  observation  is  that 
particles  reach  a  higher  peak  temperature  as  diameter  decreases  and 
the  particle  heating  time  (FWHM)  increases  as  the  particle  size 
increases.  Particles  ignited  by  the  laser  more  than  once  achieved 
more  or  less  the  same  temperatures  thereby  indicating  that  the 
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particles  did  not  burn.  No  significant  temperature  difference  was 
noted  for  particles  ignited  in  nitrogen  and  oxygen  (Fig.  28).  This 
leads  the  author  to  conclude  that  particles  with  diameters  in  the 
range  10-50  um  could  not  sustain  combustion  and  extinguished  a  few 
microseconds  after  the  ignition  pulse,  if  in  fact  they  ignited  at 
all.  The  laser  power  could  be  increased  high  enough  to  raise  the 
particle  temperature  to  nearly  2800°K.  In  spite  of  such  high  energy 
pulses,  the  particles  extinguished  without  any  significant  burning. 

4.4   ERROR  ANALYSIS 

Consider  the  sources  of  uncertainty  in  the  mass  calculation. 

m  =  K  V2/AV  , 


Differentiating, 


3V 


K  is  constant  for  all  mass  calculations  since  it  contains  only 
geometric  and  physical  constants.  Thus,  while  K  could  cause  a 
systematic  error,  it  does  not  introduce  any  uncertainty  from 
measurement  to  measurement . 

To  get  the  V  term,  V  is  measured  over  a  very  small  range  (1-2) 
at  large  voltage  (250  V).  Therefore,  the  contribution  to  the 
uncertainty  of  measurement  is  at  most  around  1%.  For  example,  in 
the  worst  case  with  a  V  »  250  and  a  voltage  range  of  2  volts, 


22  5V  =  (K/AVe)  2VSV  .  (18) 


3V 
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The  contribution  to  the  uncertainty  would  be  at  most  2%.   The 

major  source  of  uncertainty  is  the  determination  of  AV  .   Typicallv 

e    lv  -  ' 

values  of  V^  range  from  5  to  60  mv.   Since  the  precision  in 

determining  V  is  ±0.003  v,  the  effects  of  3  mv  uncertainty  can  be 

large. 

If  V  =60  mv,  uncertainty  is  5% 

If  V  =10  mv,  uncertainty  is  30% 
e  J 

The  degree  to  which  I  can  calculate  AV  is  the  key  to  the 
accuracy  of  the  measurement.  I  improved  my  certainty  of  AV  in  two 
ways.  First,  I  know  the  particle  size  to  ±2.5  urn,  and  I  know  the 
density  of  the  material.  This  gives  me  a  rough  estimate  of  the 
mass.  The  rough  knowledge  of  the  mass  prevents  gross  error  in  the 
calculation. 

The  accuracy  to  which  I  can  measure  the  temperature  by  two-color 
pyrometry  is  dependent  on  the  degree  of  accuracy  to  which  I  can 
measure  the  ratio  of  the  intensities. 

Let  the  heights  (at  a  particular  instant  of  time  with  reference 
to  the  base  linel  of  the  two  traces  be  h  and  h,  (900  and  750  nm 
respectively).   Therefore  the  ratio  of  the  intensities  is 


hl 
R=h-  (19) 


Differentiating  Eq.  (19)  and  using  only  absolute  values  we  get 


dR   dhl    <">•> 

r =  \  +  K7  (20) 


H9 


The  uncertainty  in  the  measurements  of  the  heights  could  be  off 
by  0.1  division  and  the  error  in  the  temperature  would  depend  on  the 
individual  heights  as  shown  by  Eq .  (201.  A  graph  of  ratio  versus 
temperature  is  given  in  Fig.  7. 
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5.   RECOMMENDATION  FOR  FUTURE  WORK 

Joseph  Cerv  reported  a  mass  loss  in  the  range  of  40  -  90Z  as  a 
result  of  laser  ignition.   The  author  was  however  unable  to  verify 
these  results  with  the  Millikan  chamber  for  the  following  reasons: 
i)   The  particles  (20  -  50  urn)  typically  have  a  mass  of  10   -  10 
gm.   This  results  in  a  typical  AV  in  the  range  10  -  30  mv  at 
200  V.   From  the  various  AV  obtained  it  is  very  difficult  to 
"guess"  the  AV  for  one  electron.   To  illustrate  this  point 
consider  the  "true"  AV  for  one  electron  to  be  lOmv.   If  the 
experimental  AV's  were  10,  20  and  30  mv  then  one  encounters 
little  problem  in  choosing  the  correct  AV.   Given  typical 
uncertainties  in  AV,  however,  the  experimental  AV's  would 
typically  be  7,  22  and  28  mv  thereby  making  the  task  of 
"guessing"  the  AV  for  a  single  electron  very  difficult. 
ii)   The  mass  of  smaller  particles  (lOum)  could  be  determined  by 
electron  stepping  because  AV  is  now  significantly  larger  than 
its  uncertainty,  but  these  particles  usually  did  not  produce  a 
"good"   trace   for   time-temperature  measurements.    This   is 
because  these  particles  are  too  small  and  hence  not  bright 
enough  to  be  detected  by  the  pyrometer. 

Improvements  can  be  made  in  the  present  two-color  pyrometer. 
PMT's  could  replace  the  photodiode-amplif ier  combination.  This 
would  improve  the  gain  as  well  as  eliminate  "ringing".  The  data 
collecting  system  could  be  automated  through  interfacing  with  a 
computer.   The  resolving  power  of  the  telescope  should  be  improved. 
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High  speed  cinemicrography  could  be  coupled  to  the  telescope  to  get 
the  details  of  the  "burning"  process. 


SECTION    6 
APPENDIX 
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6.1   NOTATION 

c   Speed  of  light  3  x  108  m/s 

C  Calibration  constant  of  pyrometer 

—  1  2     —  ? 
C   Planck  s  first  radiation  constant  3.74  x  10    W  cm 

C  Planck's  second  radiation  constant  1.44  cm°K 

h  Planck's  constant  6.624  x  10    J/s 

E,  Intensity  of  radiation  at  wavelength  i 

e  .  Emissivity  at  wavelength  i 

T\  Efficiency  of  detectors  at  wavelength  i 

nL  .  Efficiency  of  light  transport  at  wavelength  1 

X  Wavelength 

AA  Bandwidth 

tt  3.14159 

um  microns 

k  Boltzman's  constant 

6.2   COMPUTER  PROGRAM 

Given  a  series  of  balancing  voltages  V  ,  this  program  arranges 

the  V  's  top  to  bottom  in  an  increasing  order.  It  then  picks  up  the 

lowest  V   (topmost  number)  and  subtracts  it  from  the  rest  of  the 

numbers.   This  process  is  repeated  until  only  one  value  of  V   is 

S 

left.   The  reader  then  "chooses"  the  V   from  the  various  values 

e 

generated  by  the  program. 


PROGRAM  MILI 
C     DECLARE  VARIABLES 

INTEGER  PASS, MIN.NUM,  I, START, NEXT, DONE 
REAL  K( 100), TEMP 


74 


C    READ  IN  VECTOR 

READ(5,*)(K(I),I=1,NUM 
C    PRINT  DATA 

WRITE(7,100)NUM 
100   FORMAT ('NUMBE  OF  ELEMENTS' ',14) 

WRITE(7,200) 
200   FORMAT ('UNSORTED  VECTOR') 

DO  30,I=1,NUM 
30   WRITE  (7,*)  K(I) 

DONE=NUM-l 
C     REPEAT  FOR  NUM-1  PASSES 

DO  1,  PASS=l,DONE 

MIN=PASS 

NEXT=PASS+1 

DO  2,  I=NEXT,NUM 

IF(K(I).LT.K(MIN))MIN=I 

2  CONTINUE 
IF(MIN.NE.PASS)  THEN 
TEMP=K(PASS) 
K(PASS)=K(MIN) 
K(MIN)=TEMP 

END  IF 
1    CONTINUE 
C    PRINT  SORTED  VECTOR 

WRITE  (7,300) 
300  FORMAT  ('SORTED  VECTOR') 

DO  40,I=1,NUM 
40   WRITE(7,500)K(I) 
500  FORMAT ( IX, F J. 3) 
C    THIS  DOES  SOMETHING 

WRITE (7, 400) 
400   FORMAT ( ' EXTRA  OUTPUT ' ) 

DO  3,  START- 1, NUM-1 

DO  4,  I=START+1,NUM 

K(I)=K(I)-K(START) 
4    WRITE(7,500)K(I) 

3  CONTINUE 
STOP 
END 
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APPENDIX  6.3 
Data  sneet  for  pyrometer  components 


n 


STANDARD  BANDPASS  FILTERS 






FILTER 

TYPE           MINIMUM 

MAXIMUM 

CW/L             FWHM 

f  OF      SEE  PAGES        PEAK  T 

THICKNESS           PART  NUMBER 

INMI                IN  Ml 

CAV              6-8                      % 

N* 

IMMI 

25mm  D1A. 

50mm  DIA. 

!    H                          666.3  +21-0          10±2 

3                   3                       50 

2.0G 

7.0 

656FS10-25 

6  56  FS  10-50 

[                                660  +  3/-0             10±2 

3                   3                        50 

2.05 

7.0 

660FS10-25 

660FS 10-50 

670  +3/-0             10±2 

3                   3                       50 

2.05 

70 

670FSI0-25 

6T0FS10-50 

,    U                          670.8  +2/-0           10±2 

3                   3                       50 

2.05 

7.0 

67 1FS 10-25 

671FS10-50 

680  +3,-0             10±2 

3                   3                       50 

2.05 

7.0 

660FS10-2S 

680FS 10-50 

I                                690  +3/-0             10±2 

3                   3                       50 

2.05 

7.0 

690FS10-25 

690FS10-50 

,    Ruby  laser           694.3 +.2/-0            ]±.2 

2                   6                       60 

2.05 

3.5 

694FS01-25 

694FS01-50 

[   Ruby  laser           694.3  +.27-0           1+.2 

2                   2                       40 

2.05 

8.5 

694FS02-25 

694FS02-50 

[   Ruby  laser        _  694.3  +  .5/-0           3±.S 

2    .              2                       45 

2.05 

8.5      ... 

694FS03-25 

694FS03-50 

[  Ruby  User     _      694.3  +  1.5/-0       10±2 

3.     .'    .       3          50 

2.0S 

7-0 

694FS10-25 

694FS 10-50 

.,      700+3/-0       \    10±2 

3  .    ._V     3                       50 

2.05 

7-oV." 

700FS10-25 

700  FS  10-50 

700  +  10/-o"          40±8', 

3    .-  -.-    -.».._-.    "...  « 

2.05 

7-0 

700FS40-25 

700FS4O-50 

.   .      700 +25/- 0  '"  80±16. 

—       ...     9  ...                70 

_ 

7.0    "' 

700FS80-25 

700FS80-50 

....  .     .....  710+3/-0..'        10±2 

* T...        3 50 

2.05  . 

...  70  ... 

710FS10-25 

T10FS1G-50 

.....t  r_.  ,  ...  720+3/-70          "  10±2 

.3  ^  ,.V  *V* 50 

2.05 

V   7-o.VV. 

72OFSI0-25 

720FS1O-50 

..    ,...wL..V_73of+3/r<> Z~JQ±iZ.. 

-3___^   3  _"_ _'_'  ,50  .   ', 

2.05" 

7-0..W 

.  730  FS 10-25 

730FS10-5O 

.  .     ..-iVVVV'.^O  +3/-0^"__  10±2~" 

.jJVV.  -...3.  V  "  VV  so 

2.05 

7.0 ,.'_"   740FS10-25 

740FS10-50 

.  .  «;,  ..  — 750  +3/rO_ ,.'V.r  J0±2  VV 

- 3  V VV^VsV IV!" ."".""  50 

2.05  _ 

1 '.-..  7-o  v.v 

750FS10-25 

750FS10-50 

B   '.]t.      750 +IO1-O, "'."    40±8 - 

t        T^750+25/-o'T_~  100±20. 

T:Xi..  .„"3V,W."."'L"  "so  "* 

"-km** sVVVVVVV..™  :  _" 

2.05. 

7.0  2. 
7-0- 

750FS40-25 
750FS0O-25 

750FS40-50 
75OFS0O-50 

.-,.--_.  JJ.1_760  +31-0  V  W  10±2     ~ 

»..1"7TV.»V.  /_     .'"       « 

2.05 

M  .!/. 

760FS10-25 

760FS1O-50 

^^■jgfttfVJ.  76^5-+.2/-0.          10±2 

3     _            3        ..'.*« 

2.05 

I  7.0 .  _., 

766FS10-25 

766  FS  10-50 

al> ...._,  -  V.  .. "780  +3/r0 __'. Vl0±2.~J. 

3-V1V V  3  .V.V  VV  « 

>-lL;3V.      '    45 

2.05 
205  _ 

.......  70  VV 

7.0" 

77OFSI0-25 
780FS10-25 

T70FS10-50 
78OFS10-50 

.     .  _.."1'.790  +Z/-6 ' "io±2  "" 

Rb;j_  .  IT  _794.7 +2/-'oJ"31lO±2;.^ 
....."■  ..IVW     '  800  +3/-0"  .~I"lO±2""~" 

3" V.V«ZZV ,«.  V 

L»r  l^lu;  3 45  _ . . 

"3.  r?r3  7V.,'""« "" 

2.05' 
2.05  ' 
2.05.  . 

W- 

-.^V  7-".ViV 

790  FS  10-25 
795FS10-25 
800FS 10-25 

790  FS  3  0-50 
795FS10-50 
800FS10-50 

... ..,:■■"■  s.i.   800+10/-0       '    40±8     v 

"3;  J2Su«i  " ""  M""'" 

^^J^jC  VVVV_ !" « "  VV . 

2.05  _ 

v..  7P.zii 

800FS40-25 
800FSO0-25 

800FS40-50 
600FS00-5O 

...... . ...  TV V" 8oo.  +25/-.6IJV. ioo±2o7V 

'.  ^V_VVVV.  810,  +3/-0   .  "..V  10±2_" 

■  '■  _-  -Tlr^/ii*  820',+ 3/^*0  .ZtW.10ii.IV. 

3.JVVV  V"    JT£    45  :,. 

V»l V.~  V _..__^~l.  "«_.7 ' 

2.05  ~ 

.2.05kI 
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WVV.7  pV'XII 
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820FS 10-25 
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V.".V?oHTV 
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840FS10-25 
850FS10-25 

84  OFS 10-50 

-  -V^-—  -V—  8^°  +3XV0VVVV  10±  2 J  V 
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3=.:~V--^  3._.4.VV!'VV"5o,ZV 
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850FS40-25 

850FS40-50 
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— "V 
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FIGURE  29a 
Transmission  curve  for  tne  filter  wltn  peak  wavelengtn  750  nm. 


7P 
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FIGURE  29b 

Transmission  curve  for  tne  filter 
witn  peak  wavelengtn  900  urn. 


80 


> 

i 

\'~ 

— L_ 

=. 

j i 

- 

— 

— 

s- 

■rt 

-*- 

J— 

— 

i ,_ 

k, 

{= 

r. 

P 

§= 

r- 

j 

"n 

= 

= 

-T 

VND 

c 

■■•■ 

lEGEJJD 

r; 

- 

:=• 

'"        iSMhll 

|lfc!      ^r-X 

i 

L3T 

MS 

uu 

? 

ii 

■* 

t 

r- 

XH 

:d_ 

'c. 

- 

•= 

- 

J 

1 

i        i       i 

1      1 

:    j 

NT 

1411 

me 

i? 

l£A 

' 

j   i 

H 

3K 

1  ! 

AUC 

BD 

u 

^s 

la 

- 

-- 

i 

'0.7 

_ 

- 

-_ 

s 

-: 

; 

■1- 

I: 

£ 

= 

I 

= 

= 

= 

£ 

- 

r! 

- 

- 

- 

= 

= 

r.r 

= 

= 

ri 

= 

= 

= 

? 

~ 

:7 

\. 

~ 

.r 

.— 

= 

.= 

= 

~ 

= 

= 

= 

rrr 

.7 

?! 

J 

\ 

- 

-- 

0   0 

i;: 

:z 

3 

=.' 

rr 

§ 

~ 

= 

= 

= 

E 

- 

r 

= 

/-- 

::. 

A 

! 

r; 

r 

-r 

3 

r: 

—. 

=j 

i=E 

rrr 

= 

s 

~ 

= 

i 

= 

■-. 

= 

I: 

r-E 

= 

S 

= 

= 

= 

= 

rr: 

— 

_ 

— 

°A 

§1 

1 

= 

1 

_ 

= 

■■r-' 

N 

— 1 

— 

H 

~ 

H 

~ 

*~- 

r: 

rr 

rr 

rr 

= 

= 

= 

rr 

f 

\ 

= 

H 

=■ 

= 

= 

= 

~ 

— . 

= 

hk 

lrr 

rr 

tt 

r^ 

:^ 

\ 

0.3 

- 

=■ 

~ 

3 

i 

= 

rr 

- 

= 

— 

= 

rr 

-: 

= 

= 

-- 

\ 

= 

^.: 

:." 

r 

rr 

= 

rr. 

r_ 

~. 

i 

= 

1 

= 

= 

I 

\ 

— 

r 

= 

= 

= 

- 

E 

= 

- 

rr 

= 

rr 

= 

~ 

;-: 

\ 

I, 

K= 

= 

~ 

r 

= 

r 

rr 

— 

i_r 

/ 

? 

: 

= 

J 

'- 

™ 

\ 

O.t 

r 

rr 

V- 

~ 

• 

= 

— 

r 

_. 

/ 

- 

- 

- 

" 

-r> 

\ 



* 

— 

■5 

= 

~r 

1 

y 

5 

'-: 

i 

: 

=jg 

"|~ 

j_ 

"" 

* 

J*. 

-1 

sj 

— 

3 

~- 

-i 

-.1 

r 

\ 

J 

1 

,0,-i 

81 


S874,  S875  Series  (Visible  to  Near  IR) 


_        ,,  Outlines 

Type  No.       ,p.  q  an, 


Photosensitive 
Surface 


Size 

(mm) 


SB74  Series  (Metallic  Case) 


Effective 
Area 
(mnV) 


Package 
(mm) 


Characteristics  (®25*C) 


Spectral  Response     I  Typical  Radiant  Sensitivity  (AAV) 


Ranga 

N) 


Peak Peak  | 

Wavelength  Wave- 

(nm)      jlengih 


GaP  LED 

(560nm) 


He-Ne 
Laser 

(633nm) 


GaAs 

LED 

(930nm) 


NEP<B 

typ- 

(W/HzVi) 


typ. 

(cnvHi'-VW) 


S874-5K 


5.9x5.9 


S875  Series  (Ceramic  Case) 


To-ia 


430-1060 


S875-16R 
-33R 

0 

o 

1.1x5.9 
2.5x2.4 

5.8 
5.8 

2.7x15  I 

S50±so 

0.45 

0.2 

, 

0.3 

0.3 

6x10-" 

4x10" 

-MR 

o 

5.9x5,9 

33 

-Tixlb-430-1060 

•1010R 

• 

10.1x10.2 

98 

15x16.5) 

1x10-" 
3x10-" 

6x  10" 

3x10" 

•  Spectral  Response 


3  K :  Boroslllcate  glass  window 
R :  Resin  coating  window 
L  i  Lena  type  boroslllcate  glass  window 
i»l  nep=  Noise  Currant  (A/Hz'/i) 

Radiant  Sensitivity  at  Peak  (A/W) 
®0*  (nL.)n,,-  [Effective  SenaHlve  Area  {pm*jM 
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«  R*  =  = 10(mV) 

Dark  Currant  at  Vn  =  10m V  (A) 

nJ™XlmK9,'eqlJ"*d  ,0  ,rans'»oi  'fom  10%  to  90V.  of  the 
LED  (855nm)  and  the  load  resistance  Is  MA 
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FIGURE  31a 

Pyrometer  output  for  a  45  urn  particle  in  air. 

Laser  Power  0.30W 

Laser  Beam  Diameter  55  urn 
Time  Scale  for  this  and 

the  rest  of  the  traces  is  500  us/div 

Voltage  Scale  100  mv/div 


FIGURE  31b 
Time-Temperature  profile  for  the  above. 
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FIGURE  32a 

Pyrometer  output  for  a  42  um  particle  in  air. 

Laser  Power  0.30W 

Laser  Beam  Diameter    22  ym 
Voltage  Scale         100  mv/div. 


FIGURE  32b 
Time-temperature  profile  for  Fig.  (28a). 
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FIGURE  33a 

Pyrometer  output  for  a  30  ym  particle  in  air 

Laser  Power    0.6W 

Laser  Beam  Diameter  for  tnis  and  tne  remaining 

traces       22  ym 
Voltage  Scale  200  mv/div. 


FIGURE  33b 
Time-temperature  profile  for  tne  above. 
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FIGURE  34a 
Particle  in  figure  29a  neated  again 
Laser  Power   0.60W 

FIGURE  34b 
Time-temperature  profile  for  tne  above. 
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FIGURE  35a 

Pyrometer  output  for  a  22  um  particle  in  air. 

Laser  Power       0.40W 
Voltage  Scale     200  mv/div. 


FIGURE  35b 
Time-temperature  profile  for  tne  above. 
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FIGURE  36a 
Pyrometer  output  for  tne  previous  particle  neated  again 
Voltage  Scale         200  mv/div. 

FIGURE  36b 
Time-temperature  profile  for  Fig.  36a. 
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FIGURE  37a 

Pyrometer  output  for  a  22  urn  particle  in  air 

Laser  Power       0.30  W 
Voltage  Scale     100  mv/div. 


FIGURE  37b 
Time-temperature  profile  for  the  above. 
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FIGURE  38a 
Pyrometer  output  for  a  33  um  particle  in  air. 


Laser  Power 
Voltage  Scale 


0.30  W 

100  mV/div. 


FIGURE  38b 
Time-temperature  profile  for  the  above. 
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FIGURE  39a 

Pyrometer  output  for  a  22  um  particle  in  air. 

Laser  Power       0.25W 
Voltage  Scale     lOOmv/div. 


FIGURE  39b 
Time-temperature  profile  for  tne  above. 
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FIGURE  40a 
Pyrometer  output  for  a  25  urn  particle  in  air. 
Laser  Power       0.18W. 


FIGURE  40b 
Time-temperature  profile  for  the  above. 
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FIGURE  41a 
Pyrometer  output  for  a  "rotating"  particle  in  air. 


FIGURE  41b 
Time-temperature  profile  for  trie  above. 
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FIGURE  42a 
Pyrometer  output  for  a  30  um  particle  in  nitrogen. 
Laser  Power         0.30W. 


FIGURE  42b 
Time-temperature  profile  for  tne  above. 
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FIGURE  43a 
Previous  particle  heated  again. 


FIGURE  43b 
Time-temperature  profile  for  tne  above. 
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FIGURE  44a 
Pyrometer  output  for  a  27  urn  particle  in  nitrogen. 
Laser  Power         0.30W 


FIGURE  44b 
Time-temperature  profile  for  tne  above. 
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FIGURE  45a 
Pyrometer  output  for  a  30  um  particle  in  nitrogen. 
Laser  Power         0.30W 


FIGURE  45b 
Time-temperature  profile  for  tne  above. 
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FIGURE  46a 
Pyrometer  output  for  a  the  particle  in  Fig.  (45)  heated  again. 


FIGURE  46b 
Time-temperature  profile  for  the  above. 


129 


130 


> 


"K 


Ll 
O 

LjJ 


cc 


UJ 

_l 
u 

iH 

'HZ 
;0CUJ 

cro 
cc 


UJ 

cc 

£! 

cc 

UJ 


im 


HZ 


CEH 


UJ" 


n 

D 

z 
o 
u 
yj 
<n 

D 

cr 
u 
H 

2: 


LU 

z 

H 


1    I    §    8'    8    I    8    8    8    8 

01       co       ai        k        r-,        (i>       .■£)       in        1/1        -J- 


I    §    § 


INIA-13»«    3Unj_Ufci3dW3X 


131 


FIGURE  47a 
Pyrometer  output  for  a  25  Mm  particle  in  nitrogen. 


FIGURE  47b 
Time-temperature  profile  for  the  above. 
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FIGURE  48a 
Pyrometer  output  for  a  25  Urn  particle  in  nitrogen. 
Laser  Power         0.32W 


FIGURE  48b 
Time-temperature  profile  for  tne  above. 
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FIGURE  49a 
Particle  in  Fig.  (48)  heated  again. 


FIGURE  49b 
Time-temperature  profile  for  the  above. 
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FIGURE  50a 
Pyrometer  output  for  a  27  um  particle  in  oxygen. 
Laser  Power        0.34  W 


FIGURE  50b 
Time-temperature  profile  for  tne  above. 
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FIGURE  51a 
Pyrometer  output  for  a  35  urn  particle  in  oxygen. 
Laser  Power        0.34  W 


FIGURE  51b 
Time-temperature  profile  for  tne  above. 
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FIGURE  52a 
Pyrometer  output  for  a  30  urn  particle  in  oxygen. 
Laser  Power        0.34  W 


FIGURE  52b 
Time-temperature  profile  for  tne  above. 
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FIGURE  53a 
Pyrometer  output  for  a  17  um  particle  in  oxygen. 
Laser  Power        0.34  W 


FIGURE  53b 
Time-temperature  profile  for  tne  above. 
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FIGURE  54a 
Pyrometer  output  for  a  27  Um  particle  in  oxygen. 
Laser  Power        0.34  W 


FIGURE  54b 
Time-temperature  profile  for  tne  above. 
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FIGURE  55a 
Pyrometer  output  for  a  15  Mm  particle  in  oxygen. 
Laser  Power         0.34  W 


FIGURE  55b 
Time-temperature  profile  for  tne  above. 
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FIGURE  56a 
Pyrometer  output  for  a  30  ym  particle  in  oxygen. 
Laser  Power        0.34  W 


FIGURE  56b 
Time-temperature  profile  for  tne  above. 
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FIGURE  57a 
Pyrometer  output  for  a  30  urn  particle  in  oxygen. 
Laser  Power         0.32  W 


FIGURE  57b 
Time-temperature  profile  for  tne  above. 
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FIGURE  58 
Particle  in  Fig.  57  neated  again. 
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FIGURE  59a 
Pyrometer  output  for  a  34  um  particle  In  oxygen. 
Laser  Power         0.32W 


FIGURE  59b 
Time-temperature  profile   for   trie  above. 


170 


warn 

Hi  US!  HI 

S3H 

QMBBBK1 

inii^ 

171 


o 

LU 

J 

H 


■CD 
H 

cr 

CO 

lu; 


,_IUJ 

£UCJO 

lu*° 


cr 
cr 

UJ 
Q_ 

2: 

LlI 


LU 

I- 

H 
I 

Q. 

cr 
or 

ID 


Q 
Z 
O 

D 

UJ 

to 
o 
<r 
u 


0 

0 

0 

0 

0 

O 

0 

0 

6 

O 

O 

O 

O 

0 

0 

0 

0 

0 

a 

g 

st 

g 

Sj 

g 

i 

^ 

& 

3 

Sj 

8 

fS 

R 

ijj 

g 

Sj 

i 

$ 

I 

N 

nj 

n 

k 

fv 

(D 

<D 

IP 

tn 

* 

* 

CO 

(f) 

Oi 

0 

CNIA-13»      3Unibti3dW3X 


172 


REFERENCES 

1.  Joseph  Cerv,  M.S.  Thesis,  Department  of  Nuclear  Engineering, 
Kansas  State  University  (1985). 

2.  W.Z.  Nusselt,  Ver  Deut.  Ing.  68,  124  (1924). 

3.  D.B.  Spalding,  Fuel  30,  121  (1951). 

4.  Coffin,  K.P.  and  Brokaw,  R.S.,  Tech.  Note  3929  (1957). 

5.  Libby,  P. A.  and  Blake,  T.R.,  Comb.  Flame  36,  139  (1979). 

6.  Laurendeau,  N.M.,  Prog.  Energy  Comb.  Sci.  4,  221  (1978). 

7.  Hottel,  H.C.  and  Davis,  H. ,  Ind.  Eng.  Chem.  26,  889  (1934). 

8.  Ubhayaker,  S.R.  and  Williams,  F.A.,  Electrochem.  Soc.  123,  747 
(1976).  

9.  Frickel,  R.H.,  Shaffer,  R.E.  and  Stamatoff,  J.B.,  "Chambers  for 
the  electrodynamic  containment  of  charged  aerosol  particles", 
Chemical  Systems  Laboratory,  Aberdeen  Proving  Ground,  MD  21010 
(1970). 

10.  Arnold,  S.,  J.  Aerosol  Sci.  H),  49  (1979). 

11.  Philip,  M.A.,  Gelbard,  F.  and  Arnold,  S.,  J.  Coll.  Int. 
Sci.,  41,  507  (1983). 

12.  Hottel,  H.C.  and  Broughton,  F.P.,  Ind.  Eng.  Chem.,  4  (1932): 
166-174. 

13.  House,  K.,  B.S.  Thesis,  MIT,  Dept.  of  Chemical  Engineering 
(1978). 

14.  Dictor,  R.A.,  M.S.  Thesis,  MIT,  Dept.  of  Chemical  Engineering 
(1979). 

15.  Quinn,  T.J.,  Temperature ,  7,  294  (1983). 

16.  Ayling,  A.B.  and  Smith,  I.W.,  Comb.  Flame  _18,  173-84  (1972). 

17.  Wolfe,  W.L.  and  Zissis,  G.J.,  The  Infrared  Handbook,  2,  12 
(1978).  - 


173 


TWO-COLOR  PYROMETER  TEMPERATURE  PROFILES  FOR 
SINGLE  PARTICLE  GRAPHITE  COMBUSTION 

by 

Upendra  Brahme 


AN  ABSTRACT  OF  A  MASTER'S  THESIS 

submitted  in  partial  fulfillment  of  the 
requirement  for  the  degree 

MASTER  OF  SCIENCE 
Department  of  Physics 


Manhattan,  Kansas 
1986 


KANSAS  STATE  UNIVERSITY 


ABSTRACT 

Experimental  data  on  the  temperature-time  profile  of  laser  heated 
particles  of  AESAR  crystalline  graphite  with  diameters  in  the  range  15 
to  50  urn  in  ambient  environments  of  nitrogen,  air  and  oxygen  are 
presented.  The  modified  Hillikan  cell  was  used  to  suspend  individual 
graphite  particles.  The  temperature- time  measurements  were  made  with  a 
two-color  pyrometer.   Intensity  traces  were  recorded  to  determine 
particle  temperature  histories  and  possible  "burning"  times. 

Correlation  between  "burning"  time  and  particle  size  indicated  that 
larger  particles  "burned"  longer  and  had  lower  peak  temperatures. 
Ambient  environment  had  no  significant  effect  on  the  particle  peak 
temperatures.  Typical  peak  temperatures  ranged  from  1700  to  2200°K 
while  the  "burning"  times  ranged  from  2  to  5  ms.   Particles  ignited  by 
the  laser  more  than  once  achieved  more  or  less  the  same  temprature. 
These  results  led  the  author  to  conclude  that  particles  in  the  range  of 
15  to  50  um  could  not  sustain  combustion. 


